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Abstract 

In this work we apply point canonical transformations to solve some classes of nonautonomous 
nonlinear Schrodinger equation namely, those which possess specific cubic and quintic - time and 
space dependent - nonlinearities. In this way we generalize some procedures recently published 
which resort to an ansatz to the wavefunction and recover a time and space independent nonlinear 
equation which can be solved explicitly. The method applied here allow us to find wide localized 
(in space) soliton solutions to the nonautonomous nonlinear Schrodinger equation, which were 
not presented before. We also generalize the external potential which traps the system and the 
nonlinearities terms. 
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I. INTRODUCTION 



Investigations on the nonlinear Schrodinger equation (NLSE) [l| has increased very much 



in the last decades. The importance of such investigations are not on 
world applications of NLSE or the Gross-Pitaevskii equation (GPE) 



due the possible real 
I, but also due to the 



possibility, on the theoretical side, in increasing the class of nonlinear integrable models ^|. 

The applications of NLSE and/or GPE with spatially dependent cubic and quintic (CO) 
nonlinearities can be appreciated, for example, in pulse propagation in optical fibers |4| 
in photonic crystals and in the study of Bose-Einstein (BEC) condensates jfj|, whose 
nonlinearities are driven by means of optical interactions as well. 

Especifically, dark and bright solitons have been observed in various nonlinear physical 
phenomena. Investigation of bright and dark solitons is useful for understanding the proper- 
ties of BEC. Bright solitons are characterized by a localized maximum in the density profile 
without any phase jump across it. In the relevant experiments, this type of soliton is formed 
upon utilizing a Feshbach resonance to change the sign of the scattering length from positive 
to negative. On the other hand, dark solitons may also be considered as moving domain 
walls which separate regions of a condensate with different values of the order parameter. 
In fact, dark solitons are density dips characterized by a phase jump of the wave function at 
the position of the dip and can be generated by means of phase-engineering techniques 
Thus, there are many possibilities of managing these solitons, strongly justifying a quest for 
novel analytical solutions. 

The thecniques for managing the nonlinearities have been improved a lot and, in some 
cases, the nonlinear equations governing the system present not only space js| but time- 
dependent nonlinearities also 9|- 10). In this case one talks about the so-called nonau- 
tonomous NLSE whose localized waves solutions, in some specific cases, were found by 



Serkin and Hasegawa 



12[ resorting to a similarity transformation, which maps the nonau- 



tonomous NLSE onto a nonlinear stationary equation, whose solutions are well known. 

Similar procedure has been carried out in the recent study by Beitia, Perez-Garcia, Vek- 
slerchik and Konotop (BPVK) fl3|. Concretely, the authors relate the nonautonomous NLSE 
with cubic nolinearity (CNLSE) 



d 2 ^ o 
+ v(x,t)^ + g(x,t)\^\ 2 q, 



(1) 
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and the stationary CNLSE with constant coefficients 

/^ = -^ + G|$| 2 $, (2) 

(/i is the eigenvalue of the nonlinear equation known as chemical potential) by using the 
ansatz 

* (x, t)= p (x, t) exp (iip (x, t)) $ (X (x, t)) . (3) 

for the solution of Eq. ([1]). The variable X, which in Eq. (j2J) plays the role of space 
coordinate, is in fact a function F (£) of the specific combination £ (x,t) = j(t)x + S(t) 
of the space coordinate and the time. In such an approach, which provides analytical 
solutions for the wavefunction ty(x,t), one finds out that specific forms of the trapping 
potential v (x, t) and of the nonlinearity function g (x, t) can be explored and that there is 
an intrinsic dependence of v (x, t), g (x, t), p (x, t) and tp (x, t) on 7(t), 5(t) and F (£), that is 
one may choose conveniently j(t), 8(t) and F (£) to obtain v (x, t), g (x, t) and, consequently 
nonsingular p(x,t) and single- valued phase tp(x,t). 

More recently, in the work by Avelar, Bazeia and Cardoso (ABC) {jjj the authors followed 
the BPVK approach but extend the problem, focusing their attention on the nonautonomous 
cubic and quintic nonlinear Schrodinger equation (CQNLSE). They have also obtained an- 
alytic localized solutions of the bright or dark type (breathing, resonant, quasiperiodic and 
moving breathing solutions), depending on whether the eigenvalue p vanishes or not. 

The nonautonomous CQNLSE is obtained by adding a term h(x,t) |\1/| 4 \1/ to the right 
side of (j2J) and the stationary CQNLSE with constant coefficients is given by 

/i$ = + G 3 |$| 2 $ + G 5 1$| 4 $. (4) 

We have noticed that when the similarity transformation is applied, the trapping po- 
tential v (x,t) always has a quadratic term, namely u 2 (t)x 2 , which constitutes a harmonic 
oscillator with time- dependent frequency. The frequency is related to ^(t) and its time 
derivatives and, when d 2 5/dt 2 ^ 0, one finds a driven time-dependent harmonic oscillator, 
whose force is also time-dependent. Such kind of problem is consacred in the literature of the 
time-dependent Schrodinger equation concerning the analysis of dissipative effects in quan- 
tum fluctuations. Moreover, such systems have found real world applications in quantum 



optics 



15| and plasma physics [l6j. Among the approaches applied to solve the Schrodinger 
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equation for the time-dependent (driven) oscillator one resorts to point canonical transfor- 
mations on the coordinates and a re-scaling of the time in such a way the problem can be 
;ransformed into a Schrodinger equation for the harmonic oscillator with constant frequency 
171 ]. As a matter of fact, in some specific cases, the problem can be even reduced further 
into a Schrodinger equation for a free particle (l8|. 

One of the main goals of the present work is to apply successfully the same canonical 
point transformations for solving some time- dependent Schrodinger equations to the prob- 
lem of the nonautonomous CQNLSE mentioned above. We show that the point canonical 
transformation followed by an appropriate redefinition of the wavefunction and additional 
transformations of variables also leads to the stationary CQNLSE with constant coefficients 
(J4]), without resorting to any ansatz on the form of the wave function \1/ (x, t) in Eq. ([3]). 

The approach adopted here is straightforward and as a result of it one can see clearly 
how is the dependence of the trapping potential v (x, t) and of the inhomogeneous coef- 
ficients g(x,t) and h(x,t) with the functions 7(t), 5(t) and F (£). With this approach, 
the constraints over the functions 7(t), S(t) and F (£) in order to render the wavefunction 
\1/ (x, t) well defined are also evident. By following the mapping approach presented in 
we present the soliton solutions for Eq. (j4j) in terms of the Weierstrass elliptic function, 
such that by considering specific sets of the parameters /x, and an arbitrary constant 

of integration we not only recover the results found in [13|-[14j, but also present some new 
soliton solutions by using the set of functions j(t), S(t) and F (£) they have worked with. 

We also deal with extensions of the models by including other trapping potentials which 
are a mixing of circular functions with the time-dependent harmonic oscillator and also 
comment on a generalization of the nonautonomous time-dependent nonlinear Schrodinger 
equation by considering a non-polynomial nonlinearity. 

In the next section we present the approach to map the nonautonomous CQNLSE onto a 
stationary CQNLSE. The third section is devoted to present the solutions of fl4]) in terms of 
the Weierstrass elliptic function where some examples of wide bright and dark solitons are 
shown. In the fourth section we consider other kinds of trapping potential besides the per- 
sistent time-dependent harmonic oscillator and comment on a non-polynomial nonlinearity. 
The fifth section is left for the conclusions. 
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II. THE APPROACH 



In this section we present the approach by focusing on the nonautonomous nonlinear 
Schrodinger equation with terms of cubic and quintic order in the wavefunction, namely 



d 2 m 



dt dx 2 
whose function coefficients are written as 



+ v{x,t) V + g 3 {x,t) \ty\ 2 ^ + g 5 (x,t) |#| 4 #, 



(5) 



v{x,t) =u{t)x 2 + f l {t)x + f 2 {t)+ 1 2 {t) V[i{t)x + 8{t)\ , 



(6) 



g 3 (x,t)^G 3l (t)f[ 1 (t)x + 5(t)}, (7) 

g s (x,t) =G s h[y(t) x + 5 (t)] . (8) 

The reason for choosing the inhomogeneous coefficients in this way is going to be clarified 
below. 

Now, we implement the following coordinate transformation and time rescaling 

£ *(r) 



121 



x 



t-tn 



7 (t) 7 (t) 
dr' 



lo 7 2 if) ' 

with 7 [t (£)] = 7 (t) and S [r (t)} — 8 (t) . Then, one can recast (jSJ) as 



(9) 
(10) 



-I 2 



d 2 ^ 
W 



v(Z,T)y + G 3 ^f(Om 2 V + G 5 h (0 |^| 4 



(11) 



where 7 r = dj/dr , 5 r = d<S/dr , \P (£, r) = \P (as (£, r) , t (r)) and 



(12) 



7 7 

From the last two equations one can appreciate why we have chosen the specific dependence 
of V , / and h on £ = 7 (t) 2 + 5 (£). Eq. (iTTjl looks like a nonautonomous NLSE on £ and 
r, except for the first derivative term in the variable £. In order to eliminate that term, we 
redefine \I/ (£, r) as 

^^vWe-^^^r) , (13) 



5 



where a (£, r) = j=£ 2 + (d T — -^j | — a (r) , with a (r) an arbitrary function for the moment. 
By substituting (ITS]) in (TTTj) one gets 

^ = -^ + ^(e,r)^ + G 3 /(OI^N + G 5 / i (OI^|V, (14) 



where 



PttlT)=(B(T) _^iia! + f 7lW _^^ 



4 y 7 4 y^' 2^7 



One can see why the factors involving 7 (t) are present in the expressions of v (x, t) and 
#3 (x, i). Moreover, one can also appreciate the contribution of the redefinition ()13|) . specif- 
ically the contribution of the phase a (£, r) to the redefinition of the trapping potential 

Now, one can choose 7 (r) and 5 (r) such that 

and, by coming back to the original variables (x,£), the functions w, fi and /2 appear as 

U {t) ~ 4 7 2 ' 2? ' /2() ""4y~^' (16) 

revealing the intrinsic connection between the frequency u (t) , the force fi (t) and the 
functions 7 (t) and Note that a(t) is an arbitrary function that could be chosen as 

da/dt = (<5f/27) 2 if /2 (t) were not present in ([6]). Thus, one removes the explicit time- 
dependency of dHJ), that is 

< = -|| + v (0 ^ + g 3 / (0 H 2 V + c 5 (0 H 4 V>- (17) 

and the wavefunction f)13p is written as 



* (x, t) = VTWe-'"^') V (£ <) , I" (*)) , (18) 

where a (x, t) = ^ x 2 + ^x — a(t). We recall that the phase a (x, t) was introduced through 
the redefinition (TT31) and it contributes to the redefinition of the trapping potential, which 
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could be eliminated thanks to the presence of the time- dependent driven harmonic oscillator 
terms (the three first terms v (x,t)). 

For stationary solutions in the variable r, that is r) = (£) exp(— iEr), we have 

^ = (V (£) - E) + G 3 f (0 |0| 2 + G 5 A (0 |0| 4 • (19) 

Since we still have a nonlinear equation with inhomogeneous nonlinearities, we are go- 
ing to make further transformations in order to reach a nonlinear second-order differential 
equation with constant couplings. For that we redefine £ as a function of another variable 
C, that is £ = F (£), which lead us to the differential equation 



d 2 6 F, 



CC<1<P -FUV[C(0]-E) 



d( 2 F ( d( 

+G 3 f 2 f [C (C)] |0| 2 + G 5 F 2 C h K (C)] |0| 4 , (20) 
where = dF/dC, and = (£ (£)). By redefining the field as 

0(0 = ^^(0, (21) 

we reach Eq. (jSJ) in terms of ( 
where 



-/x$ + G 3 |$| 2 $ + G 5 |$| 4 $, (22) 



2 

£1 3F 1 1 

v = -fI(v[U0]-e) + -^-^-, /[£(C)] = =3, Me(C)] = =4- (23) 



In terms of the variables £ and r we have 



V ^ = (w) 2 -^ F ' 2 + E ' ^) = F ' 3 ' M£) = ^\ (24) 

where F' = dF/dl; and the field (l2ip becomes 

0(O = ^=$(C(O)- (25) 

Finally, by returning to the original space-time coordinates (x,t), the wavefunction can 
be obtained from (TTHj) and f[2"5"j) . that is 



* (x, t) = exp [-i V t)} $ [F (£ (*, t))] , (26) 



with 77 (2, t) = j^x 2 + j^x — a(i) + E J Q * cfa'7 2 (t'), where a(t) is an arbitrary function. 

Thus, we have shown, by means of transformation of variables, how the nonautonomous 
CQNLSE, Eq. can be mapped onto a nonlinear one-dimensional second-order differential 
equation with cubic and quintic nonlinearities Eq. 0221) . More than that, we have shown 
explicitly how the part of the trapping potential V(£) and the nonlinearities functions /(£) 
and h(£) are related to the transformation function = F (£) (see eqs. I\Z£ 



III. THE MAPPING ONTO THE WEIERSTRASS p-FUNCTION 

Here we consider the cases in which $ is a real function on £ (122 p . One can check that 
the solutions for eq. ( I22p are also solutions of the equation 

2 

d( J " *~ ~ ' 2 ~ ' 3 
where e is a real arbitrary constant. Moreover, we shown how to map such kind of equations 



M ) =e-/i$ 2 + %$ 4 + ^$ 6 . (27) 



19| onto the nonlinear Weierstrass differential equation 



= ^p 3 -92p-gs, (28) 



where g 2 and g% are the Weierstrass invariants and the discriminant is A = g\ — 27 'g\ . The 
values of the invariants and of the discriminant determine how the Weierstrass function are 
written in terms of the double-periodic Jacobi elliptic functions for A / 
solutions for ( 12 7p can be included in one of the categories listed below 
1 ) For e 7^ one has the following solution for $ (C) 



20|, and the 



$(C) ^ P (C. feg3 )H-,/3 - <29) 

where g 2 = |/i 2 — 2Gse and g% = ^/j 3 — |G*3/ie — |G 5 e 2 which leads to the following 
discriminant A = -±f (6eG^ + 36e 2 G 2 + 36G 3 G 5 e/i - 3^ 2 G 2 - 16/i 3 G 5 ) . 

2) For e = and \x 7^ one sees that g 2 = |/i 2 , g% = J^/i 3 and A = 0, then the solution 

is 



where (5 = ±y ^ + |/i G 5 . 

3) For e = and /i = one gets #2 = #3 = A = (in this case p (£) = (~ 2 ) and 

2 Gfi ^-A2\-V2 



*(C) 



G 3 



(31) 



3G 3 2 

All the examples with real world applications considered in jl3| and jl^j may be repro- 
duced and fit in one of the these categories. We present next three additional examples 
which were not considered by them. 



A. Examples 

Because the dark and bright solitons are of primary importance for developing concrete 
applications of BEC, many technique are developed for manipulating and controlling the 
soliton's parameters and induce changes in their shapes which would be useful for applica- 
tions. One possibility is to vary the atomic scattering length by means of external magnetic 
fields, i.e. by using Feshbach resonances. People has been asked to themselves on how to 
modulate the width and amplitude of a soliton in a controllable manner. It has been demon- 
strated that the variation of the scattering length provides a powerful tool for controlling 
the generation of bright and dark soliton trains ; especially, this technique can be used to 
modulate the bright soliton into very high local matter densities in both harmonic trap 



potential 



2l| and repulsive (inverted) harmonic trap potential [22 1 



We show by means of some examples that the width and the amplitude of the solitons can 
also be manipulated theoretically by focusing on the same form of the functions F (£) , 7 (t) 
and 5 (t) considered in Jl^| but with some different values of /i, G3, G5 and e. Explicitly, we 



take F (f ) = ^sErfi (^^J (this is the imaginary error function |20|), 7 (t) =- y 1+3cos 2 (2f) 
and S (t) = 0. We remark that this choice of F (£), when substituted in V{^) (HM)) . implies 
into an extra contribution for the time-dependent frequency u(t) of the driven harmonic 
oscillator, or in other words, the trapping potential v (x, t) is a sole time-dependent driven 
harmonic oscillator for fi = 0. 
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1. Example 1: 



In real world applications the nonlinearities and the confinement potential can be changed 
independently using an optical trap and magnetic-field- induced Feshbach resonance. In order 
to modulate the width of the dark soliton in a controllable manner, we gradually increase 
the self-interaction and simultaneously turn off the trapping potential. 

Here, there is no way to turn off the ubiquitous trapping harmonic potential, but we can 
modulate the width of a dark soliton through a suitable choice of the chemical potential and 
the coefficient of the quartic self-interactions. Here, as in real world applications, they are 
related to each other. 

First we notice that Eq. (1221) can be seen as the static equation of motion d 2 §/d( 2 = 
dU/d$, where U($) = (l/2)(-yU$ 2 + + ^ $ 6 + c). If the parameters of the potential 
are chosen such that it presents only two minima, say at ±1, and such that the chemical 
potential varies from positive to negative values (we notice that d 2 U/d§\ =Q = —/J,), the 
minimum energy solution $(0, which connects the minima of {/($) at ( — > ±oo, is a kink 
which is deformable into a two-kink (double kink) as fi varies from positive to negative 
values. 

Such $ 6 polynomial potentials have been used to study phase transitions which comes 



with domain wall splitting and the appearance of a wet phase in some ferroelectric 23| and 
paramagnetic J^J materials. The wetting transition can also take place in deconfinement 
phase transition of SIK3) Yang-Mills theory (25 1, in supersymmetric QCD 2(| and in thick 
brane- world scenario J27 1, (28 1. 

We choose fi = 2a 2 — 1 , G3 = 2 (a 2 — 2) , G5 = 3 and e = c = a 2 . The so- 
lution $(C) ; f° r the chosen coefficients of {/(<&), falls into the first category of solu- 
tions shown in ( I29p with g 2 = | (1 + a 2 ) 2 and g 3 = — (1 + a 2 ) 3 and p (£ ', g 2 , #3) = 
I (1 + a 2 ) + (1 + a 2 ) csch 2 (Vl + a 2 £) , such that we have 



_ tanh (Vl + a 2 () , 
$ (C) = ^ =y= . (32) 

'sech 2 (VTT^C) + a 2 



The profiles of |$ (C)| 2 f° r three different values of \i are shown in Fig. ([Q), from which one 
can see the increasing of the width of the soliton as fi becomes close to the critical value 

/!=-!. 



10 



By using Eq. (12tj|) we find the wavefuncction 



7 2 x 2 



Qb 2 



sinh ( y/l 



X 2 V3 
2G 



exp [—i 7] (x, t)] x 
1/3 ^n j 



1 + a 2 cosh 2 (v/TT^ Erfi (^l) ) 
which is a wide breathing dark soliton for a 2 < 1/2 . 



I<i>|- 





- v=fj/ 




i \ \ 0.8 






1 \ 1 °' 6 






1 \ K.4 










' ,V=-0.9998 




(33) 



FIG. 1: Left: Profiles of |$ (£)| from Eq.®- Right: |* (re, t)| for the wide breathing bright soliton 
with A = 0.001, b = 10, a = 0.01, 



2. Example 2: 

Bright soliton can also be modulate into a desired width and amplitude in a controllable 
manner by changing the scattering length and the trapping potential. 

As in example 1, the soliton width can be modulated by a convenient choice of the 
parameters. Here we take values of [i less than the critical value /x = — 1 of the previous 
example and U ($) no longer positive, but with two symmetrically disposed global minima 
and one local minima at $ = 0. We look for solutions such that $(£ — > ±oo) — > 0. For that 
we take e = c + k = 0, where k is a constant of integration. From relations 18.12.1 to 18.12.3 
of 20| one can see that p (£, 4/i 2 /3, 8// 3 /27) = /3)(1 + 3 csch 2 -\/j/i[C) and by choosing 

11 



2/3/G 3 = -A 2 < (that is, G 5 < 3G 2 /16 and G 3 < 0, such that y/ — ~/G~ 3 is real, we 
find that 



A 2 cosh (2^7^) + 1 



*(0 = |i + ifp(C,y2,^)-^ ) 



We show profiles of $ (() for three different values of A 2 < 1 and \fj,\ = — G3 = 4 in Fig. 
(J2J) . In the same figure we show also the amplitude |\P (2, t)| 2 for the wide breathing bright 
soliton. 

Particularly, a very thin bright soliton can be obtained when G5 = 0. We set \i = — 1 
(the critical of /i), G3 = — 1 and e = , such that #2 = § , #3 = _ Jp, A = and $ (£) 
belongs to the second category of solutions listed in the previous section, namely 

-i\ -V 2 

2 ' 2 \ r 3 

One can verify that p (Cs^S^) = | + csch 2 £ and, consequently $ (£) = -\/2sech£. 
By using the Eq. (126]) we arrive at the breathing bright soliton solution 

*{x,t)=e*' 6 yffiexp(-?^ expH77(x,t)]sech^-^|^Erfi^^j . (36) 

where a(t) = J (7* (E - £)) dt 

3. Example 3: 

This is an interesting example because we can find three different analytic solutions with 
the same set of parameters fi = 5 , G3 = 10 and G5 = —3 . Among those solutions we find 



periodic solutions which were not presented in 



In the case e = and using ( 130]) we find two, very similiar to each other, periodic solutions: 



$1 (C) = vl ° and $ 2 (C) 



(5 - V5) + 2\/Esm 2 (VEC) J {5 + V5) - 2A/5sin 2 (v^C) 

(37) 

However, if e = 5 ( 4v ^ 5 ) we g n( j ^he kink-like configuration 
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FIG. 2: Left: Profiles of |$(C)| 2 from Eq.(j5I]). Right: \^ (x,t)\ 2 for the wide breathing bright 
soliton with A = 0.001, b = 10 



$(0 = 




which provides a breathing dark soliton solution for \1/ (x,t). 

IV. GENERALIZATIONS 

In this section we discuss some generalizations of the nonautonomous nonlinear 
Schrodinger equation. Particularly, we focus on two generalizations. 

A. Generalization of the trapping potential 

As it was shown in equations f lM|) . part of the potential is determined by a convenient 
choice of F(^), such that F' (£) does not have zeroes. Moreover, one can see that V (£) 
has, for fi = E = 0, a strucuture similar to the one- dimensional supersymmetric quantum 
mechanics potential, with F" /2F' playing the role of supersymmetric superpotential. Then 
it would be interesting to analyse the influence that several trapping potentials in susy 
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quantum mechanics could have while keeping the same nonlinearities on the nonautonomous 
nonlinear Schrodinger equation. 

We have analysed some of those possibilities and have found that the trigonometric Scarf 
and trigonometric Rosen-Morse potentials which have in their spectrum only bound-states 
(here we are using the same nomenclature for susy quantum mechanics potentials as that 



of |29|) are good candidates for the entrapment of breathing bright solitons in a spatially 
periodic trapping potential. 

For the sake of brevity we just present the first case here, that is we take 

F' (0 = G~ 1/3 sec 2 ^ «) , (39) 

where A, a > and -tt/2 < a£ < tt/2. Then V (f) = A (A - a) sec 2 (a£) - A 2 - 
fiG~ 2/S sec 4A / a (af) and g 5 (x,t) = G 5 G~ 4/3 sec 8A / a «). 

It is always convenient to look for 7 (t) and S (t) such that ^ (x, t) is well defined. We 
choose 

7 (t) = 1 + l + 7l sin (t) + 72 sin (V2t) 2 , 5{t) = 0. (40) 

30| for more details on this function 7 (t). The time-dependent frequency and 
driven force are given by 0161) . 

In order to get explicitly brightlike solution, let us consider /i = 0, G3 = 2, G5 = —3 and 
e = 0. The solution is given by 

^ = 2 1 /y/ 2 exp (x,t)]cos A /" K) (41) 

v 1 + F (0 2 

whose modulus squared is shown in figure ([3]) for A = a = 1. 

A breathing dark solution is get by setting /x = 3, G3 = 6, G5 = —3 and e = 0. It is given 

by 

^ = 6 l/6 7 l/2 exp [_j g (X; F (0 CQS A/ a (<) 

/ 1 + F(0 2 



We refer to 



B. Generalization of the nonlinearities 

We have shown in equations (124|) how the nonlinearities functions g% (x, t) and (x, t) (see 
Eqs. ©-©) must be intrinsically related to each other and given in terms of the functions 
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FIG. 3: |^ (x,t)\ 2 for the bright soliton from Eq. (|4ip with 71 = 0,72 = 0.1. The periodic structure 
is due to the periodicity of the circular trapping potential. 

7 (t) and F' (£), such that the mapping of Eq. (jSJ) onto fl22|) could be realized. From a close 
inspection on those relations we have found that the approach used here can be used to 
map a nonautonomous Schrodinger equation with non-polynomial nonlinearity onto a static 
sine-Gordon like equation. In the first paper in Ref. 30| it was proposed a polynomial 



nonlinearity up to Nth order in t), with fuctional coefficients related to each order such 
that the nonautonomous NLSE could be mapped onto a stationary NLSE with constants 
coefficients. Here we present inverse transformations to obtain a nonautonomous NLSE 
with non-polynomial nonlinearity from a stationary NLSE with non-polynomial nonlinearity. 
Particularly, we consider the sine-Gordon equation which is satisfied by $(C)> namely 



q^ = ^ lbnQ]= ±^ b ^» ( 43) 

where b is positive constant. By using the transformation ( = F (£) and the redefinition of 

(44) 

Now, by considering that 0(£) = ip (£, t) e lEr , we find that ^ (£, r) satisfies the nonau- 
tonomous NLSE 



$ (F (0) = VF' (0 <p (0 we find that 

\ '\ 00 ( 1 \ n 

pi \n+2i2n-l i2n+l 
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where V ® = (g,)' E. 

Now, the inverse of transformation in equations f[9"j)- (lTUj) . namely 



£ = 7 (t) x + 6 (t) and t - r = / dt'Y (f ) 

'o 



(46) 



lead us to 



7 2 <9r 7 3 <9x 7 2 <9x 2 

oo /_-i\n 

^ 2n + l ! v ; m ^ 



n=0 



(47) 



where ^ = (x,£) , r (£)). The first-derivative in x can be eliminated by the following 
redefinition of the wavefunction ip(x,t) = e %a ^ x,t ^ (x,t) /\/^y (£), where ct(x,t) = ^x 2 + 
— a (t). Then we finally obtain the nonautonomous NLSE with a generalized nonlinearity 



.3^ <9 2 # 



dt 



dx 2 



oo fiY" /p / \ n+2 



5H 



5W*(x,t) 



(48) 



where 



v (x, t) 



7«7 - 27 2 2 ( <5«7-27A 



4 7 2 



ar + 



2 7 2 



x 



and 



(49) 



H 



2 



d 2 



+ v(z,t) 



+^7 S * V [£(*,*)] cos 
In this specific generalization the soliton amplitude is periodic and given by 



i*(x,or 



47 (t) arccos ( 1/ sn [F(^)/\/b\ 



(50) 



(51) 



where sn (z\u) is the snoidal Jacobi elliptic function with elliptic parameter < v < 1. For 
v = 1, sn(z\v) = tanh(z) and we have a breathing dark soliton. 
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V. CONCLUSIONS 



In this work we have shown that an approach used to tackle the Schrodinger equation with 
explicit time-dependent potential parameters can be extend to the case of non-autonomous 
NLSE. As a consequence we were able to reproduce results originally obtained through the 
use of an Ansatz together with similarity transformation. Besides, this procedure lead us 
to novel solutions of that problem and, particularly some configurations of wide soliton 
solutions. 

Furthermore, we extended the approach in order to deal with systems with generalized 
nonlinearities and trapping potentials which are a mixing of the (driven) time-dependent 
harmonic oscillator and circular functions. Extensions of this work for the cases with two 
or more fields are under development. 
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